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Abstract 
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I. INTRODUCTION 



One of the most convincing manifestations of the difference between the BCS scenario 
and superconductivity in the cuprates is the pseudogap, or the depletion of a single particle 
spectral weight around the Fermi level |J. This is observed mainly in the underdoped 
cuprates where the pseudogap opens in the normal state as the temperature T decreases 
below the crossover temperature T* and extends over a wide range of T. 

Due to the complex nature of cuprate systems, there are a number of theoretical ex- 
planations for the pseudogap behavior. One of them is based on the model of a nearly 
antiferromagnetic Fermi liquid 0. Another possible explanation relates the pseudogap to 
spin- and/or charge- density waves ||. A third direction, which we take in this paper, argues 
that precursor superconducting fluctuations may be responsible for the pseudogap phenom- 
ena. Indeed an incoherent pair tunneling experiment |4]] proposed recently may allow one to 
answer whether the superconducting fluctuations are really responsible for the pseudogap 
behavior. Furthermore, one cannot exclude the possibility that the pseudogap is the result 
of a combination of various mechanisms, e.g., both spin and superconducting fluctuations. 

The precursor superconducting fluctuations have recently been extensively studied using 
different approaches. In most cases, the attractive 2D or 3D Hubbard model was considered. 
In particular this model has been studied, both analytically |§-|7| and numerically [f^TT 



in the conserving T-matrix approximation that is "$ derivable" in the sense of Baym |T2| 



The non "$ derivable" T-matrix approximation was considered in [I3|. In this approach, 
the pseudogap is related to the resonant pair scattering of correlated electrons above T c . For 
the d-wave pairing, the pseudogap was also studied in fl4| (for a review, see [|L5|0 and Monte 



Carlo simulations for the 2D attractive Hubbard model were performed in [16[ 



It is known, however, that while the T-matrix approximation provides an adequate de- 
scription of 3D systems at all temperatures, including the superconducting state with a 
long-range order, it fails (see, for example, 0) to describe the Berezinskii-Kosterlitz-Thouless 
(BKT) transition into the state with an algebraic order, which is only possible in 2D sys- 
tems. This is why, in most of the papers cited above, the T-matrix approximation was used 
to study either 3D systems 0,^,^,0 or 2D systems above T c ||^, |I^|I3J in order to avoid 
the BKT transition, even though it is generally accepted that 2D models are more relevant 
for the description of cuprates JI7|] . 

Of course, the superconducting transition itself is not of the BKT type, because even 
a weak interplanar coupling produces a transition in the d = 3 XY universality class, suf- 
ficiently close to the transition temperature. Outside the transition region, however, the 
low-energy physics is governed by vortex fluctuations |18| , and one can expect the 2D model 
to be especially relevant for the description of the pseudogap phase. This was confirmed for 
the quasi-2D model |19| (see also |2( 



Regarding the pseudogap, it is sufficient to consider the case where T > T c . However, one 
definitely needs an approach different from the T-matrix if one wants to study the 2D theory 
for the entire temperature range and wants to connect the pseudogap to the superconducting 
gap. An alternative approach overcoming the above difficulty was proposed in |21^p3|. For a 
2D system, one should rewrite the complex order field $(x) in terms of its modulus p(x) and 
its phase 9(x) as $(x) = p(x) exp[i9(x)], which was originally suggested by Witten in the 
context of 2D quantum field theory [p4fl . It is impossible to obtain $ = ($(x)) 7^ at finite 
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T because this would correspond to the formation of symmetry breaking homogeneous long- 
range order, which is forbidden by the Coleman — Mermin — Wagner — Hohenberg (CMWH) 
theorem p3 |. However, it is possible to obtain p = (p(x)) ^ with $ = p(exp[i6(x)}) = at 



the same time because of random fluctuations of the phase 9{x) (i.e., because of transverse 
fluctuations of the order field originating in the modulus conservation principle |^6|). We 
stress that p ^ does not imply any long-range superconducting order (which is destroyed 
by phase fluctuations) and, therefore, does not contradict the abovementioned theorem. 

For the simple model studied in [pl] , p2| , there are three regions in the 2D phase diagram. 
The first one is the superconducting (here, BKT) phase with p ^ at T < T BKT , where T B kt 
is the BKT transition temperature, which plays the role of T c in pure 2D superconducting 
systems. In this region, there is an algebraic order or a power law decay of the ($*$) 
correlations. The second region corresponds to the so-called pseudogap phase (Tbkt < T < 
T p ), where T p is the temperature at which p is supposed to become zero. In this phase, p 
is still non-zero, but the above correlations decay exponentially. The third is the normal 
(Fermi-liquid) phase at T > T p , where p = 0. Note that $ and all the symmetry violating 
correlators like ($(s)$(0)) vanish everywhere. 

The proposed description of the phase fluctuations and the BKT transition is very similar 
to that given by Emery and Kivelson [27]. However, the field p(x) does not appear the 
phenomenological approach of [p7[] , while in the present microscopic approach it occurs 
naturally. We also mention here the application of similar ideas to the 3D case ||28|| , where 
instead of the 2D temperature Tbkt one has the phase transition temperature in the 3D 
XY-model, T c xy . 

The main quantity of interest in the present paper is the one-fermion Green's function 
and the associated spectral function A(u>, k) = — (l/ii)lmG(uj + i0, k). The second quantity, 
being proportional to the intensity of the angle-resolved photoemission spectrum (ARPES) 
29fl , encodes information about the pseudogap and quasiparticles. Following the approach 



of Refs. |2~T|-p3"| the Green's function for the charged (physical) fermions is given by the 
convolution (in momentum space) of the propagator for neutral fermions (which has a gap 
p 7^ 0) and the Fourier transform of the phase correlator (exp(ir 3 9(x)/2) exp(— ir 3 9(0)/2)) . 

Thus, the approximation employed here assumes the absence of coupling between spin 
and charge degrees of freedom; this can be taken into account at the next stage of approx- 
imation. We demonstrate that the quasiparticle spectral function broadens considerably 
when passing from the superconducting to the normal state, as observed experimentally 
|2T| . More importantly, the phase fluctuations result in a non- Fermi liquid behavior of the 



system both below and above T B kt- 

We note that the effect of classical phase fluctuations of the order field on the spectral 
properties of underdoped cuprates has also been analyzed by Franz and Millis ||30|| . Being 
experimentally motivated, they could show that the corresponding photoemission and tun- 
neling data are well accounted for by a simple model where <i-wave charge excitations are 
coupled to supercurrent fluctuations. 

A brief overview of the paper is as follows: In Sec. O, we present the modulus-phase 
formalism for the fermion Green's function and explain why it is so important to use this 
formalism for the description of 2D models. In Sec. [TIT], we obtain and discuss the Green's 
function of phase fluctuations both below and above Tbkt- This expression is then used in 
Sec. [TV] to derive the temperature and retarded fermion Green's functions. We show that 
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this Green's function exhibits a non-Fermi liquid behavior. In Sec. [Vj, we obtain an analytic 
expression for the spectral density of the fermion Green's function and discuss this result in 
detail. The density of states (DOS) is considered in Sec. [VI]. Appendix [A| contains technical 
details on the calculation of the long-distance asymptotic behavior of the phase correlator. 
Appendix [B] contains the derivation of an alternative representation for the fermion Green's 
function which is useful in calculating the spectral density. The integrals for the DOS are 
given in Appendix 



II. THE MODULUS PHASE REPRESENTATION FOR THE FERMION 

GREEN'S FUNCTION 



Our starting point is a continuum version of the two-dimensional attractive Hubbard 
model defined by the Hamiltonian density |2l|-|23 



(2.1) 



where x = r, r denotes the space and imaginary time variables, ip a {x) is a fermion field with 
the spin o =t? I, Tn is the effective fermion mass, fi is the chemical potential, and V is an 
effective local attraction constant; we take h = ks = 1- The model with the Hamiltonian 
density (|2.1|) is equivalent to the model with an auxiliary BCS-like pairing field, which can 
be written as 



n = ^H 



r 3 



2m 



r_<$>*(x) 



ty(x) + 



|$(x) 
V 



(2.2) 



in terms of Nambu variables 



y(x) 



^t(x) 
tfj\(x) 



* t (a;) = (^(a;) ^(x) ) 



(2.3) 



where r± = (t± ± ir<i)l% T3 are the Pauli matrices and $(x) = V^{x)r^{x) = Vipiip] is 
the complex order field. 

We consider the full fermion Green's function in the Matsubara finite temperature for- 
malism 



G(x) = (*(x)*t(o)) 



(2.4) 



For the 3D case of the BCS theory, the frequency-momentum representation for ( 2.4 ) in the 
mean field approximation is known to be [3"I| 



G(iuj n , k) 



iuj + r 3 ^(k) - r + $ - r_$* 
a;2 + e(k) + |$| 2 



(2.5) 



where us n = (2n + l)irT is the odd (fermion) Matsubara frequency, £(k) is the dispersion law 
of electrons evaluated from the chemical potential /1, and $ = ($(x)) is the complex order 
parameter. 



4 



A problem arises when one tries to apply Eq. ( |2.5|) directly to 2D systems, since it has 
been proved (see p5| ) that nonzero <3> values are forbidden. Nevertheless, one can assume 
that the modulus of the order parameter p = |$| has a nonzero value, while its phase 9{x) 
defined by 

$(x) = p{x) exp [iO(x)\ (2.6) 

is a random quantity. To be consistent with (|2.6|) one should also introduce the spin-charge 
variables for the Nambu spinors 

= exp[ir 3 0(x)/2]T(x), #t( x ) = T+(z) exp[-fr 3 0(x)/2], (2.7) 

where T is the neutral fermion field operator. The strategy of treating charge and spin 
(neutral) degrees of freedom as independent seems to be quite useful, and at the same time 
a very general feature of 2D systems p4T , |32"|| . 

Applying fl2.7|) , we thus split the Green's function ( j2.4|) into spin and charge parts 



Gap{x) = E GM*)({^* y2 )M^ {0)/2 )rp), (2.8) 

a', /3' 

where 

G*p{x) = (T a (x)Tj(0)) (2.9) 

is the Green's function for neutral fermions. Introducing the projectors P± = |(J ± r 3 ) we 
obtain 

e ir 3 0/2 = p +e ^/2 + p_ e -W/2 j e — 3^/2 = p_ e tf/3 + p +e -W/2 j (2.10) 

so that (|2.8| ) can be rewritten as 

G(ar) = ]T i 3 a^(a;) J P/3(exp[^(x)/2]exp[-z^(0)/2]), (2.11) 

where a = (3 and a = —f3 correspond to the diagonal and non-diagonal parts of the Green's 
function, respectively. 

For the frequency-momentum representation of ( 2.11 ) we have 

00 f d 2 p 

G(iu n ,k)=T 22 TT~v2 E p <xG( iuJ m , p)PpD ct j 3 (iu; n - iu m , k - p), (2.12) 

m=-oo J y^ 71 ) a,P=± 



where 



rl/T r 

G(iu m ,k) = J dr J d 2 r exp[iu m r — ikr]£/(r, r) (2.13) 



and 



rl/T f 

D aP (in n ,q) = dr J d 2 r exp (itt n T - iqr) (exp[ia0(r, r)/2] exp[-i/?0(O)/2]) (2.14) 



is the correlator of phase fluctuations with even (boson) frequencies Q n = 2imT. 

There is a good reason to believe (see [p2]Q that for T close to T B kt, the fluctuations 
of the order parameter modulus p (the so-called longitudinal fluctuations, which in fact 
correspond to carrier density fluctuations and undoubtedly must be taken into account in 
the very underdoped region) |33| are irrelevant and one can safely use the Green's function 
( 12.13] ) of the neutral fermions in the mean- field approximation (compare with ( |2.5|) ) 



re i \ lu nI + r 3 £(k) - np , , 

g(zu n ,k) = 2 T C 2(\ \ _i_ 2 2 - 15 



Here £(k) = k /2m — p with k being a 2D vector and p = (p(x)). Note that in [f21| , p2jl , p{x) 
was treated only in the mean-field approximation, which means that fluctuations in both 
p{x) and 6{x) were neglected, and therefore a second-order phase transition was obtained 
at T p . However, as stressed in Introduction, experimentally the formation of the pseudogap 
phase does not display any sharp transition and the temperature T* observed in various 
experiments is to be considered as a characteristic energy scale, rather than as a temperature 



where the pseudogap is reduced to zero f3~4fl . We believe that taking the p(x) fluctuations 



into account may resolve the discrepancy between the experimental behavior of T* and the 
temperature T p introduced in the theory. 



III. THE CORRELATION FUNCTION FOR THE PHASE FLUCTUATIONS 



As stated above, we expect the phase fluctuations to be responsible for the difference 
between properties of the charged and neutral fermions defined above. The latter are de- 
scribed by the Green's function (|2.15| ), which coincides with the BCS Green's function (|2.5|) 
only under the assumption that the phase 9 of the order parameter $ = pexp(i8) is a con- 
stant and can be chosen to vanish. This is not the case for the 2D model, where there is a 
decay of the phase correlations and the Green's functions of charged and neutral fermions 
are nontrivially related via Eq. ( [2.12| ). To establish their relationship, one must know the 
correlator for the phase fluctuations. Its calculation is quite straightforward for T < Tbkt, 
while for T > Tbkt one can apply the results of the BKT transition theory . 



A. The correlator for T < T B kt 

In the superconducting phase, the free vortex excitations are absent and the exponential 
correlator is easily expressed in terms of the Green's function 



D e (x) = (0(x)0(0)) 
here, as above, x = r, r) via the Gaussian functional integral 

•VT. r.n rl 



exp 



/ T>6(x) exp 

1 rl/T rl/T r r 

2 Jo Jo dl ~ 2 ■/ ri / rf2r2/ ( ri > ri ) jDe ( ri ~ T2 ' ri ~ r 2K( r 2,r 2 ) 



(3.1) 



(3.2) 
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with the source 



I{x 1 ) = - 
The Green's function 



.^(n - r)5( ri - r) + z^n^n), (a, fi = ±). 



Dg L (x) = -J(p,T,p)V 2 r -K(p,T,p)(d T 



(3.3) 



(3.4) 



for this model was found in |22 . Note that the superfluid stiffness J and compressibility K 
are here the functions of p, T and p, and also that the Green's function (|3.4j) includes only 
the lowest derivatives of the phase 9. The higher terms are also present in the expansion, 
but we neglect them. In the simplest case where J(/i, T, p) ~ n/, the density of carriers, and 
K(fjL,T,p) ~ const f22f . 



Substituting (|3.4|) into ( |3.2| ), we obtain 



D a p{x) = exp 



OO 

-- y 



d 2 q 1 — a/? cos(qr — O n r) 



(2tt) s 



Jg 2 + KQl 



(3.5) 



It is easy to see that for zero frequency Q n = 0, the integral in Eq. (|3.5|) is divergent at 
q = unless a = (3, and therefore only two terms survive in the sum over a, f3 in Eq. ( p,12|) , 
namely 



p-g{iuj n , k)p_ + p+g(iu n , k)p + 



iu; n J + r 3 ^(k) 
'a;2+e 2 (k)+p 2 ' 



(3.6) 



It is important that the terms like P±G(iu n , k)P T , which are proportional to T\ and thus 
violate the gauge symmetry, do not contribute to Eq. ( |2.12| ) due to vanishing of the cor- 
responding -D+_ and D |_ correlators standing after them. This explicitly demonstrates 

that the non-diagonal part of the 2D Green's function is vanishes at all finite temperatures. 
Thus, making use of the Gor'kov equations for the calculation of its diagonal part and the 
gap function is questionable. For nonzero correlators, we have 



D(x) = D A 



D (x) = exp 



T 



E 



qdqdip 1 — cos(gr cos <p) cos Q n T 



(27T) 



Jq 2 + KQ 2 n 



(3.7) 



In what follows, we consider in detail only the static case r = 0. The restriction to this 
case is one of the few main assumptions we use throughout the paper. 

The summation over n and the integration over ip in ( |3.7|) can be readily done yielding 
the following exponent of ( |3.7| ) 



16nVJK 
where we introduced the scale 



J dqe- q/A [l - J (gr)]tanh 



qrg 
4 



(3.8) 



TV K' 



(3.9) 
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which is a function of the variables used (in the simplest case, r ~ \fn^/T). In ( p.8|) , we 
introduced the cutoff A by means of the exponential function. This cutoff represents the 
maximal possible momentum in the theory, the Brillouin momentum. 

One can derive from (|3.8| ) (see Appendix [A]) the following asymptotic expressions 

D(0,r) ~ < Vr ° 7 T (3.10) 
I (tt)" 8 * 7 ' r»A- 1 »r . 

This long-distance behavior governs the physics of ^-fluctuations that we intend to study in 
what follows. 

We now discuss the meaning of the value r . Using again the phase stiffness J(T = 0) 



and compressibility K from |22| we readily obtain that r = 2yep/m/T, which is the single- 
particle thermal de Broglie wavelength (ep — irnf/m is the Fermi energy). Then, assuming 
that T ~ T B kt and taking T B kt — £f/% |EI|,E3, we can estimate 



16 W2 

r ~ — = = — — , 3.11 

y/eFm kp 

where kp is the Fermi momentum. The value of kp for cuprates is less than the Brillouin 
momentum A, which is why the first case in (|3. 10 ) seems to be more relevant. 



There is another way to estimate r^. we can use the value 2A/T C , and hence, 

r ~ V2n— ^o, (3.12) 

where = Vp/ijrA) is the BCS coherence length. This shows that r has the meaning 
of a coherence length, which appears to be rather natural since the minimal size the phase 
coherence region should be of the order of £o- Since the coherence length in cuprates is larger 
than the lattice spacing A -1 , we again obtain that the first case in ( |3.10| ) applies. Therefore, 
for T < Tbkt and for static fluctuations, we have that 

T 



D(t)= (^) ^ (3-13) 



where r = ldf-^/epm. 



B. The correlator for T > T BKT 

For T > Tbkt, the expression for static correlator ( 3.1 3|) can be generalized using the 
well-known results of the BKT transition theory ||5 , 3t| , 



D(r) 



in 



J 



exp 



(3.14) 



where 
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M?) = C exp J J" T T (3.15) 

V 1 — 1 BKT 

is the BKT coherence length and C is a constant whose value is discussed later. One can 
consider Eq. (|3.14 ) as a general representation for D(r) for both T > Tbkt and T < Tbkt 



if the coherence length £+(T) is considered to be infinite for T < Tbkt- The pre-factor in 
Eq. (|3.14|) is related to the longitudinal (spin-wave) phase fluctuations, while the exponent 
is responsible for the transverse (vortex) excitations, which are present only above Tbkt- 
The pre-factor appears to be important for a non-Fermi liquid behavior discussed in what 
follows. Note, however, that the longitudinal phase fluctuations can be suppressed by the 
Coulomb interaction [3(J that is not included in the present simple model. One further 
comment is that while the approximation used to study the vortex fluctuations in is 



good for T well above Tbkt ; the form of the correlator D is appropriate for T close to Tbkt- 
The constant C can be estimated from the condition that £+(T) cannot be much less 

than the parameter r which is a natural cutoff in the theory and we thus take C = r /4 in 

our numerical calculations C = r /4. In any case, for T>Tbkt, where ( |3.15|) is valid, the 

value £+(T) is large and not very sensitive to the initial value of C. 

There also exists a dynamical generalization of ( |3.14j ) proposed from phenomenological 



backgrounds in [37 



T 

D(t, r) = exp(-Tt) (—) 8nJ exp (--1—) . (3.16) 



Note that t is the real time and 7 is the decay constant, and theherefore ( |3.16| ) is the retarded 
Green's function. We hope to consider the more general case of dynamical phase fluctuations 
( |3.16| ) elsewehere. 



C. The Fourier transform of D(r) 



For the Fourier transform ( |2.14j ) of ( |3.14j ), we have 

-l/T 



D(iQ n , q) = j drj cfVexp (ift n T - iqr) (r/r )- T ^ J exp(-r/^ + (T)) 



Z7T— — r, 



T/8nJ 




1 

drr 1 - T/87Tj J (qr) exp(-r/^+(T)). 



T Jo 

The integral in Q3.17 ) can be calculated (see, for example, [pql ) with the result 



(3.17) 



D(iQ n ,q) 



6 n , 2rcr 2 {1 - a) r(2a) 



2F1 [a, -a + -; 1 



2' 'qi + (l/Z + )\ 



(3.18) 



The hypergeometric function F(a, b; c; z) can be well approximated by a constant since 
it is slowly varying at all values of q. As this constant, we can take the value of the 
hypergeometric function at q = 00. Thus, 



o n ,o 



D(in n ,q) = ^A[q* + (l/Z + y\ 



2i-a 



(3.19) 
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where 



A = 



47iT(a) 
r(l-a) 



2 N 2(a-l) 

r , 



a = 1 — 



T 



16ttJ 



(3.20) 



It should be stressed that for T > Tbkt, the parameter a quickly deviates from unity 
as €p decreases; in other words, the underdoped region has to reveal highly non-standard 
properties in comparison with the overdoped one. 

Note that for C,^ 1 = (T < Tbkt), Eq. Q3.19|) is an exact Fourier transform of the 
correlator ( [3.13|) . 

One should take into account that even for T < Tbkt, the propagator ( ft.l9| ) does not 
have the canonical behavior ~ 1/g 2 , which is typical, for example, for the Bogolyubov mode 
in dimensions d > 2. In 2D, the modes with a propagator ~ 1/q 2 would lead to severe 
infrared singularities [^5|; to avoid them, these modes transform into softer ones (~ l/(? 2c \ 
a < 1). 



Finally, substituting ( |3.6| ) and (|3.19| ) in Q2.12| ), we obtain 



G(iu n , k) = -A J 



d 2 q iuj n + r 3 ^(q) 



(2tt)2 u>1 + ^(q) + P 2 [(k - q) 2 + ' 



(3.21) 



The coincidence of the Matsubara frequency in the left- and right-hand sides of Eq. (p. 21 ) 
is evidently related to the static approximation used in this paper. As we see in the next 



sections, the Green's function (|3.21| ), spectral density, and the density of states can be 
evaluated exactly. 



IV. THE DERIVATION OF THE FERMION GREEN'S FUNCTION 

The fermion Green's function can be calculated analytically if we split the fermion part 
of ( pD as 



where 



At 



+ 



An 



iuj + r 3 £(k) 

UJI + ^(k) + p 2 " e(k ) + iJ u 2 + p 2 • e(k) _ i 1^2 + p 2 



(4.1) 



A 1 = -|r 3 



UJr. 



Using the representations 



1 °° 

= =F« / dsexp [±is(a ± ib)] 

a±ib J 



(4.2) 



(4.3) 



T(a) 



a—l —ct 



dtt a ~ L e 



(4.4) 
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and taking ( |4.1| ) into account, we can rewrite (|3.21 
%A 



as 



G(iu n , k) 
W 1 



A\ exp 



/ ds / dtt a - 1 e*+ t-W^+p 2 x 
io Jo 

A 2 exp 



9 



js — ius — (k — q) t 

2m 



q 

-is h ius — (k — q) 2 t 

2m 



(4.5) 



Note that the special form of the integral representation ( |4.3| ) (compare with representation 
( |4.4| )) guarantees that the Gaussian integral over q is well-defined independently of the sign 
°f £(o) = q 2 /2m — /i. Now the Gaussian integration over momenta q in ( |4.5|) can be done 
explicitly: 



47iT(a) jo jo 



4vrr(a 



A 2 



A, 



St 



i — is /2m 



exp z 



2m t — is /2m 



IflS 



t + is /2m 



exp 



2m i + is /2m 



+ 2/iS 



• (4-6) 



Changing the variables as s -> 2ms and further as t — > st, we can integrate over s with 
the result 



G(iu> n ,k) 



imA r°° 



2tt Jo 



dtt 



a-l 



Atit-i) 



a-l 



£+ t(t + 2m Jul + p 2 (t-i)- itk 2 + 2imfx(t - i) 



A 2 (t + i) 



a-l 



i+t{t + i) + 2mJul + P 2 (t + + itk2 ~ 2imfi(t + i) 



(4.7) 



In the general case where ^ l ^ 0, the denominator of (|4.7|) is quadratic in t and some 
further transformations are needed. Replacing t — > — iu and expanding the quadratic poly- 
nomial in the denominator, we have 



G(iu n , k) 
where 



2a 



27T 



U [(u + u 1 )(u + u 2 )} a Jo 



-ioo A 2 u a - 1 (u + 1) " 1 
du-. 



[{u + u^u + u 2 )Y 



(4.8) 



/ fc 2 £ 2 _|_ 1 , \ 

ui = mi 2 + 1 2 ^2 + + P 2 + ^ ) 



with 



D 



u 2 = m£ 2 



k 2 e + + 1 

2m£ 2 



fc 2 e + 1 

2m£i 



+ i^l + P 2 ) + —pi^-^ul + p 2 ) 
m^_, 



(4.9) 



(4.10) 



and 



1H = MilV^n + P 



< + P 2 



(4.11) 
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We can verify from ( |4.9|) that Reuj > for /i < 0, and therefore, we can rotate the 
integration contour to the real axis: 



G(iuj n , k) 



Am!; 



2a 



A^ju + iy- 1 rr—2 

du- - — + (V^n + P 



2n [Jo [(u + u\)(u + u 2 ) 



+ (4.12) 



The integral representation ( |4.12| ) can then be analytically continued to p. > 0. The change 
of the variable z = u/{u + 1) allows Eq. ( |4.12| ) to be expressed in terms of Appell's function 



Fi(a,/3,/3',7;x,7/) 



r(7) 



1 z^Hl - zY'- a 



-a-l 



r(a)r(7 -a) Jo (1 - zxf{l - zyf 



■dz, 



(4.13) 



and hence, 



G(iw n , k) 



AmC,_ 



2a 



2na 



Ai ( u x - 1 u 2 - 1 

-b\ a, a, a; a + 1; 



Ml«2j V «1 «2 

+ ( V /^T7^- V /^T7) 



(4.14) 



For T < Tbkt; the BKT coherence length is infinite = 0), which means {u\ — 1)/mi = 
1 in the first argument of the Appell's function. This allows us to apply the reduction formula 




F, (a, (3, (3', 7; x, 1) = r(7)r(7 « g) ^ ( ^ ^ ^ 

T(7 - a)r(7 - /?') 



(4.15) 



and express the result via the hypergeometric function 

/ r, \ a-l 



G(iu n ,k) = -T 2 (a) 



At 



A 2 



-(/i-W^ + P 2 )] c 



iFt a, a; 1; 



k 2 /2m 



1-(p + *M + p 2 )Y 



- 2 F X a, a; 1; 



k 2 /2m 



(4.16) 



where we inserted the value of A from (|3.20|) . 

This completes our derivation of the temperature fermion Green's function. 



A. The retarded fermion Green's function 

To obtain the spectral density, we need to obtain the retarded real-time Green's function 
from the temperature Green function by means of analytical continuation iuj n — > uo + iO, 



and where yu 2 +p 2 -^ iy/uj 2 — p 2 . This results in the following rules (compare with ( |4~2l 

(PI), (CT)) 





A ^ A *='^-73h? )' (417) 
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U\ — > vi = mE s+ ^ p. — \Ju) z - p 2 + VV 



with 



u 2 -> v 2 = mil f ^Hg 1 - !->■ - \/-> 2 " P 1 ~ Vvj, (4. 18) 



D ^ v= \ i i^-^-^n + 4 ( " +v/ ^ T7) (4l9) 

and 

Vi = Vi(y[tf - p 2 -> -^U 2 -p 2 ). (4.20) 

For the retarded Green's function we thus have 

G(u;, k) = — — < — 1\ I a, a, a; a + 1; 



27ra [ (t> 1^2)° V f 1 t>2 



+ ( ^2_ p 2^_^2_ p2) ^ (4.21) 

It is easy to see that 



Vl v 2 = -2m£\{ii + y/^-p 2 ). (4.22) 

We now discuss the condition under which the imaginary part of G(u + iO, k) is nonva- 
nishing. 

For |co>| < p, we can see that V\ = v^, v 2 = v 2 , and therefore G(u,k) is real and 
ImG(u; + iO, k) = 0. The case where \u\ > p is more complicated. It follows from the 
Appell's function transformation property f4(| 

F 1 (a, (3, (3', 7; x, y) = (1 - z)-^ f a, 7 - - /?', /?', 7; -^-r, ^) ■ (4.23) 

V x — 1 1 — x J 

that for real x and y, the function i<\ becomes complex if x > 1 or/and y > 1. This implies 
that G(u,\s) has an imaginary part if V\ < or/and t>2 < 0. Looking at the expressions 
Q4.18Q for Vi and t>2, we can see that V\ is always positive, while v 2 may be negative. This 
means that G(u,la) has a nonvanishing imaginary part if v±v 2 < 0. Using ( J4.22 ), the 
condition for the existence of a nonzero imaginary part of G*(u;,k) can then be written as 
p + y/u' 2 - p 2 > 0. 



B. The branch cut structure of G(uj,\s) and a non- Fermi liquid behavior 



We now consider the retarded fermion Green's function ( |4.16|) for T < Tbkt- Applying 
the analytic continuation rules from the previous subsection to Eq. ( 4.16| ), we obtain 



a-l 



G{u,k) 



-T 2 (a) 



mr£ 
Ao 



A x 



k 2 /2m 



Fi a, a; 1; 



k 2 /2m 



p + yjuj 2 - p 2 ) 



-(p - V^7)} a2Fl r a] 1; a* - V^7), 



(4.24) 
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Near the quasiparticle peaks where u ~ ±E(\t), the arguments of the hypergeometric 



function in ( 4.24 ) are close to 1. One can consider, for instance, the first hypergeometric 
function, then 



k 2 /2m 



(4.25) 



Using the relation between the hypergeometric functions 1 40| 



T(c)T(c — a — b) 
2 F 1 (a, h- c; z) = ) ' K jriF^a, b; a + b + 1 - c; 1 - z) 



+ 



T{c-a)T{c-by 
T( c)T(a + b-c) 

mm 



[l-z 



\C—a—b 



2 F 1 (c-a,c-b;c+l-a-b;l-z), (4.26) 



we obtain that near z\ 



a-l 



G{u,k)~ -r 2 



[a 



r(l-2a) r(2a-l) 



mrn 



[-(^ + v/^- P ^)]«[r 2 (i 



a) 



r 2 i 



a 



\2a-l 



(4.27) 



It can be seen that the expression for the Green's function obtained is evidently a nonstan- 
dard one: besides containing a branch cut, it clearly displays its non-pole character. The 
latter in its turn corresponds to a non-Fermi liquid behavior of the system as a whole. It 
must be stressed that non- Fermi liquid peculiarities are tightly related to the charge (i.e., 
observable) fermions only, because the Green's function ( |2.15|) of neutral fermions has a 
typical (pole type) BCS form. In addition, it follows from fl4.27p that new properties appear 
as a consequence of the ^-particle presence (leading to a ^ 1), and because the parameter a 
is a function of T (see ( |3.2U| )), the non- Fermi liquid behavior is developed with temperature 
increase and is preserved until p vanishes. 

It is interesting that in Anderson's theory Kl], it was postulated that the Fermi liquid 
theory is broken down in the normal state as a result of strong correlations. Here, we 
started from the Fermi liquid theory and found that it is broken down due to strong phase 
fluctuations. As suggested in ETJ , the non- Fermi liquid behavior may lead to the suppression 



of the coherent tunneling between layers, which in turn confines carriers in the layers and 
leads to the strong phase fluctuations. In contrast to |41|], however, our model predicts the 
restoration of the Fermi liquid behavior as T decreases, since a — > 1 as T — ► (see the 
discussion in Sec. |V G| item 4.). 

The T = limit can also be obtained as follows. Strictly speaking, one cannot estimate 
the value of ro in the limit as T — - > in ( J4.24 ) via Eq. ( |3.11| ), because the substitution 
of Tbkt — € f/8 in (|3.9|) is not valid in this case. However, this is not essential because 
T/8nJ — > 0, so that the correlator ( |3.13|) , D(r) — ► 1, which evidently means the formation 



of a long-range order in the system. Furthermore, the value of a in Q3.20D goes to 1 as 
T — > 0, and the hypergeometric function in ( 4.24j) reduces to the geometrical series, 



iFAlA-l-z) 



1 



(4.28) 



Therefore, inserting (|4.28 ) in ( |4.24j ), we obtain the standard BCS expression 
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G u (u,k) 



u> + e(k) 
u 2 -£ 2 (k) - p 2 



(4.29) 



for the diagonal component Gn(uj,k) of the Nambu-Gor'kov Green's function G(u,k). 

Evidently Eq. fl4.29|) results in the standard BCS spectral density |31| with two (^-function 
peaks 



A{u,k) = - 



1 + 



E{k) 



6{u - E(k)) + - 



E{k) 



5{u + E{k)), 



(4.30) 



where E(k) = y£ 2 (k) + p 2 . To recover the nondiagonal components of G, one has to restore 
the correlators _D_ + (r) and _D + _(r) that were omitted in Sec. |III A. 



V. THE SPECTRAL DENSITY OF THE FERMION GREEN'S FUNCTION 



As is well known, |j3Tj| , the spectral features of any system are entirely controlled by its 
spectral density 



A(u, k) = — ImG u (^ + i0, k) 



7T 



(5.1) 



which, for example, for cuprates is measured in ARPES experiments (see |[29||). This function 
defines the spectrum anisotropy, the presence of a gap, the DOS, etc. In what follows, we 
calculate A(u, k) for the Green's function obtained above. 

A. Analytical expression for the spectral density 



For Vi > and v 2 < 0, the retarded fermion Green's function ( [4.21] ) can be rewritten 
(see Appendix [FJ) as 



G( W ,k)=-^[^ 

1 r(i-a) 



■l) a T(a)T(l-a) ( v 2 (l - Vl )\ 

-2-Fi a, a; 1; — - | + 



Ml-^ 2 )] c 



vAl - v 2 ) 



] — \T~7~, \ F ± f 1 ' a, 1 - «; 2 - a; — , — 

\v 2 \L{2 — a) V v 2 u 2 



u 2 -p 2 )\. (5.2) 



Then, according to (|5.1|) the spectral density for the Green's function (|5.2|) has the form 



. . Am(T sin(7ra) „. 2 2 . 
A(u, k) = — -sgnu; 0(uT - p 2 ^ 



x 2 F 1 a, a; 1; 



2tt 2 
v 2 (l - Vi 



(A 



i ii- 



r(a)r(i-a) 

[vi(l -v 2 )] a 



2—^ J 9(ji + ^ 2 - p 2 ) - ( ^ 2 - p 2 - ~^ 2 -p 2 ) 



Vx{l - v 2 ) 

Using the quadratic transformation for the hypergeometric function 



(5.3) 



2 F 1 (a, b; a - b+ 1; z) = (1 - z)-\F x ( ~, -b + 1 + a - b; 



'1-z) 



(5.4) 
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the expression ( |3.20|) for A, Eqs. ( |4.18j ), and (|4.19| ) we finally obtain 



a-l 



A(u,k) 




sgnu; 6{uj 2 — p 2 ) x 



,F X I 1; -4 2mV " ' ^" ^ ] (9(/x + V^ 2 - P 2 ) 



(5.5) 



where the chemical potential p can be, in principle, determined from the equation that fixes 

the carrier density |[22|| . Here, however, we assume that the carrier density is sufficiently 

high and p = sp- 
in the BCS theory, A(u, k) given by Eq. (|4.30| ) consists of two pieces that are the spectral 

weights of adding and removing a fermion from the system respectively. Note that our 

splitting of A(u, k) is different since each term in (|5.5|) corresponds to both the addition and 

the removal of a fermion. 

In the next subsections, we verify the sum rule for (|5.5|) , plot it for different temperatures, 

and discuss the results. 



B. The sum rule for the spectral density 



It is well known that for the exact Green's function G?(u;,k), the spectral function (f>A 
must satisfy the sum rule 



dojA(u,k) = l. (5.6) 
The Green's function calculated in ( 4.21) is, of course, approximate. This is related to the 



use of the long-distance asymptotic behavior ( |3.10|) of the phase correlator ( |3.7|) . This means 
that its Fourier transform ( 3.19Q is, strictly speaking, valid for small k only, while we have 



integrated our expressions to the infinity. Another approximation that we have made was 
the restriction to the static phase fluctuations. Thus, it is important to check whether the 
sum rule ( |5.6[) is satisfied with sufficient accuracy. 

It is remarkable that for ( |5.5|) , the sum rule ( |5.6|) can be tested analytically with the help 
of the techniques used in calculating N(u) in Appendix We obtain 

duA(uj,k) = -^t- (5.7) 
) 1 [2 — a) 

The numerical value of the integral at the temperatures of interest can be estimated as 
follows. Taking the phase stiffness J = 2/7tT B kt at T = T B kt, the value a from ( p.2U|) is 
given by 

T~Tbkt (5.8) 

J* J- BKT 
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for T close to Tbkt- In particular, a(T = Tbkt) = 31/32 gives the following estimate for 
the right-hand side of (|5.7f) , T(a)/r(2 — a) ~ 1.037. This shows that for T ~ T B kt, the 
spectral density ( pT5| ) is reasonably good at the temperatures of interest. 

The parameter a can however differ strongly from unity at T > Tbkt and in the under- 
doped regime. 



C. Results for the spectral density 

The plots of the spectral density A(u, k) given by (|5.5| ) at temperatures below and above 
Tbkt are presented in Figs. 1-3. To draw these plots, we used the value of a from Eq. ( |5.8| ) 
and the mean-field value of p obtained from the corresponding equation in [^l],^. From 
these figures and our analytical expressions, we can infer the following results: 

1. For T < Tbkt (the case presented on Fig. 1), there are two highly pronounced quasi- 
particle peaks at uj = ±E(k). They are simply related to the contribution of zeros of 
V (see Eq. §T|)) to A(u,k). 



We also observe two peaks at uj = ±p when k ^ kp (for k = k F , the value E(k F ) = p, 
so that the two sets of peaks coincide). One can check that the divergence at these 
points is weaker than at the former peaks at uj = ±E(k). In fact, these peaks are the 
result of the static and nointeracting approximation for the phase fluctuations used 
here. They are essential to satisfy the sum rule (|5.6| ). 

If the dynamical fluctuations are taken into account, it is clear that the "external" 
frequency uj in v4(u;,k) is different from the "internal" frequency in ^4 1; ^4 2 (see the 
discussion after Eq. ( |3.21| ) and compare it with Eq. ( |2.12| )). We believe that this 



additional summation over the "internal" frequency (which is present if the dynamical 
fluctuations are considered) would considerably smear these peaks, moving the excess 
of the spectral weight inside the gap. This assumption is supported by the results 
of |37| (see item 3 below). The same effect can also be reached when the interaction 



between the charge and spin degrees of the freedom is taken into account | 30| . Note 
also that the full cancellation of these peaks takes place in the T = case given by 
Eq. flO|). 

3. For uj < \p\, we have A(cu, k) = and a gap exists at all T (including T > T B kt)- This 
result is also a consequence of the static approximation used above. The dynamical 
fluctuations should fill the empty region resulting in the pseudogap formation in the 
normal state. Indeed, a filling of the gap was obtained in a related calculation |37 



where the correlator (exp(i8(r, t) exp(— i8(0))) (which differs from ( |3.16| ) only by the 



factor 1/2 muliplying the phase), which includes the dynamical phase fluctuations, was 
used in the numerical calculation of the self-energy of fermions and in the subsequent 
extraction of the spectral function from the fermion Green's function. 

In the approximation used in the present paper, the spin and charge degrees of freedom 



are decoupled (see Eq. (|2.8|) ). However, this coupling can be included at the next stage 
of approximation and also leads to a pseudogap filling. Indeed, using the special form 
of the scattering rate proposed in |j39| , it was obtained in |3(| that A(uj,h) ^ even 
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for uo < \p\. On the other hand, as stated above, there also are indications |37] that 
a filling of the gap can be obtained by considering the dynamical phase fluctuations 
only. At present, it is not clear which of these gap filling mechanisms plays the main 
role; this is the subject of our current investigations. 

4. The main peaks at u = ±E(k) have a finite temperature dependent width which is, 
of course, related to the spin-wave (longitudinal) phase fluctuations. As T — ► 0, the 



width goes to zero, but this limit cannot be correctly derived from (|5.5|) because this 



is an ordinary function, while the BCS spectral density (|4.30|) is a distribution. The 



correct limit can, however, be obtained for the integral of A(u, k) (see Sec. [VI|, where 



the density of states is discussed). This sharpening of the peaks with decreasing T in 



the superconducting state was experimentally observed [29[] and represents a striking 



difference from the BCS "pile-up" ( |4.30[ ) which are present for all T < T c 



It was pointed out in |30| that the broadening of the spectral function caused by these 
fluctuations can be greater than the experimental data permits. This leads |J0] to 
the conclusion that the spin-wave phase fluctuations are probably suppressed by the 
Coulomb interaction. 

5. For T > Tbkt (see Figs. 2,3), one can see that the quasiparticle peaks at uj ~ ±E(k) 
are less pronounced as the temperature increases. Indeed, the value of A(uj, k) at 
uj = ±E(k) is, in contrast to the case where T < Tbkt, already finite. This is 
caused by the fact that T> ^ since £ + is already finite due to the influence of the 
vortex fluctuations. As the temperature is increases further, £ + decreases, so that the 
quasiparticle peaks disappear (compare Figs. 2 and 3). This behavior qualitatively 
reproduces the ARPES studies of the cuprates for the anti-node direction [^J (see 
also [^]) which show that the quasiparticle spectral function broadens dramatically 
when passing from the superconducting to normal state. 

6. It is important to stress that due to a very smooth dependence of ^ on T (see Eq. 
( |3.15|) ) as the temperature varies from T < T B kt to T > T B kt, there is no sharp 
transition at the point T = Tbkt- There is a smooth evolution of the superconducting 
(excitation) gap A sc = p into the gap A PG , which also is equal to p and in fact can 
be called a pseudogap because the system is not superconducting at T > Tbkt- This 
qualitatively fits the experiment |2^,[54],[|!| and appears to be completely different from 
the BCS theory ]3~T |, where the gap vanishes at T = T c . As was already mentioned, 



the gap obtained at T > Tbkt takes place in the static approximation only and begins 
to be filled after dynamical fluctuations are taken into account (see, for example, [fj7[| ). 

7. Again for T > T B kt, one has A(uj, k) = when |u;| < p, and we expect the gapped 
region to be filled by the dynamical phase fluctuations [ 3~7| . We predict, however, an 
essential difference between the filling of the gap at T > T B kt and T < T B kt- Indeed, 
due to the presence of the vortices above Tbkt, the value of the decay constant 7 in Eq. 



( |3.16| ) should be much larger than for T < Tbkt- This and a nonzero value of above 



Tbkt may explain the break at T = T c in the scattering rate T 1 introduced in P9fl . In 
general, it is interesting to establish a correspondence between the phenomenological 
parameters, Ti and T introduced in [p^] and the vortex parameters £ + and 7 used 
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here. Note, however, that this correspondence cannot be simple because of the non- 
pole character of the Green's function derived here. 

As mentioned above (see item 3), the filling of the gap due to dynamical phase fluctu- 
ations is not the only possible mechanism for filling and the presence of vortices above 
Tbkt can be taken into account via coupling the spin and charge degrees of freedom 
[ p0|1 . It could also be that both these mechanisms are physically equivalent, since they 
relate the gap filling to the presence of vortices in the system. 

Since we used the mean-field dependence p(T), it is clear that the distance between 
the quasiparticle peaks (which is approximately equal to 2p) diminishes as T increases. 
This process of the pseudogap closing is accompanied by the destruction of the quasi- 
particle peaks. It is evident also that for p = 0, the normal Fermi liquid behavior 
is immediately restored because J(p = 0) = [pTjp2| . Recall, however, that the de- 



scription proposed here cannot be applied when p is rather small, because, as already 
mentioned, the fluctuations of p(x) have to be also taken into account in this region. 



VI. THE DENSITY OF STATES 



The density of states can be found from the formula 

r d 2 k r w k 2 



(6.1) 



where Nq = m/2-7r is the density of 2D states in the normal state (W is the bandwidth). 
This integral can be calculated analytically (see Appendix |C|), and which gives 



N{u) 



(A 



Hi 




a-l 



sgno; 9(ui 2 — p 2 ) x 



1-a 



1-a 



+ W-p 



(v^-p 2 ^-^ 2 -^ 2 )}- 



uj 2 — p 2 



6(p + Jlu 2 - p 2 



Again for T = and large /i>/), Eq. ( |6.2|) reduces to the BCS result |31 | 



N(u) = N 



^Juj 2 - p 2 



(6.2) 



(6.3) 



The plots for DOS fl6.2|) are presented in Fig. 4 (T < Tbkt) and Figs. 5,6 for T > Tbkt, 
respectively. If one does not pay attention to a small difference in the curves shown in these 
figures, it is well seen that qualitatively the form of DOS does not differ from standard BCS 
curves. Moreover, similarly to the spectral function, the DOS in the static approximation 
has a gap both above and below Tbkt and does not reveal any change when the temperature 
crosses the phase transition point. This confirms once more the crossover character of the 
latter, although, as was already pointed out, a 2D system is superconducting below T B kt 
only. According to generally accepted views, the existence of an empty gap above the 
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critical temperature is impossible. The reasons for its persistence were discussed in the 
previous section. Recall only that dynamical fluctuations or fluctuations of the modulus p 
undoubtedly result in the gap filling above T B kt- One must also take the dependence of the 
decay constant 7 into account (see Q3.16Q ), which for T > T B kt can be essentially bigger 
than in the region T < Tbkt due to the presence of the vortices. 

From the physical point of view, the filling of a gap (transforming it into a pseudogap) 
above Tbkt (or T c in quasi-2D case) has to continue up to T* (or T p if there is a point 
where p = 0). However, taking p-fluctuations into account (i.e., p(x) —* p + Ap) will cause 
the appearance of the self-energy, in addition to p 2 , in the denominator of the mean field 
Green's function ( |2.15| ); it is proportional to the quantity (Ap(x)Ap(O)) whose contribution 
could persist at all T. In this case, the beginning of the pseudogap opening will be defined 
by the experimental technique sensitivity of the spectral function or DOS measurement. 



VII. CONCLUSION 

To summarize, we have derived analytic expressions for the fermion Green's function, 
its spectral density, and the density of states in the modulus-phase representation for the 
simplest 2D attractive Hubbard model with the s-wave nonretarded attractive interaction. 

While there is still no generally accepted microscopic theory of HTSC compounds and 
their basic features (including the pairing mechanism), it seems that this approach, although 
in a sense phenomenological, is of great interest since it enables one to propose a reasonable 
interpretation for the pseudogap phenomena related to the vortex fluctuations. The results 
presented here are entirely analytic, which allows a deeper understanding than in the case 
of a numerical investigation. In particular, the analytic investigation of the Green's function 
structure revealed that the phase fluctuations lead to a non-Fermi liquid behavior below and 
above Tbkt- 

Evidently, there are a number of important open questions. The main question is whether 
the pseudogap is related to some kind of superconducting (in our case, phase) fluctuations. 
Hopefully, the experiment proposed in Q may answer this question. It seems plausible 
from the theoretical point of view that superconducting fluctuations should contribute to 
the pseudogap (see, however, [0). Nevertheless, one cannot exclude the possibility that 
the superconducting contribution may be neither the only nor the main contribution. 

Another open question is which approach allows one to obtain the pseudogap from the 
attractive Hubbard model. The schemes used in []3D| and in our paper are very different from 



those of ]T3[ . In particular, our approach allowed us to establish a direct relationship between 



the superconducting fluctuations and the non-Fermi liquid behavior in a very natural and 
transparent way. Also, it relates the pseudogap to the "soup" of fluctuating vortices (see also 
|30|,43l)) while [13] emphasises the existence of metastable pairs above T c . It is possible that 



both these pictures capture some physics, but in the different regions of the temperatures. 
When T is high and close to T*, the value of p is small, so that p-fluctuations or metastable 
pairs dominate. Then as the temperature approaches Tbkt, the values of p and the phase 
stiffness J are growing bigger, so that the vortex excitations dominate and p-fluctuations 
become less important. We stress once more that the vortex excitations cannot be adequately 
described within T-matrix approximation M. 
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Recently the last part of this picture was supported experimentally |44| by the measure- 
ments of the screening and dissipation of a high-frequency electromagnetic field in bismuth- 
cuprate films. These measurements provide evidence for a phase-fluctuation driven transition 
from the superconducting to normal state. 

Finally, there remains the problem of a more complete treatment of the pseudogap in 
the modulus-phase variables. In particular, the effects of dynamical phase fluctuations and 
the fluctuations of the order field modulus must be considered. The latter are especially 
important for the ci-wave superconductor since the modulus can be arbitrary small in the 
nodal directions. In this case it will be important again to check the complete structure 
of the Green's function, especially its non-pole structure. Another important question that 
has to be addressed is which factor is more important for the gap filling, the spin-charge 
coupling proposed in |3(| or dynamical phase fluctuations which, as was shown in [[J7[], also 
result in filling. 
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APPENDIX A: THE ASYMPTOTIC OF THE PHASE CORRELATOR 



To calculate the integral in Eq. fl3.8|), we first write it as 

1= / e- q/A \l - J (gr)]cothag = 
J o 

1 r°° / 
- / dte~ t/Aa cotht- 
a jo V 



1 - Jn -t 



+ 1 - r^e-^[l-J (rt)]. 
a Jo t 



(A.1) 



The following formulas are used when calculating J: 



„ * e -*(coU. t -I)=lnf + I-*(l + J 

r°° dt 



r°° at 

\ T e- pt [l - J (ct)] = In 
Jo t 



p + y/p 2 + C 2 
2^ ' 



(A.2) 



Hence, we obtain 



7 = 1 

^0 

4 



1 + y/1 + (Ar) 2 Aro / 2 

111 A 1" —A W 1 + 1 

Ar 4 V Ar 

, r Ar / 2 
r 4 V Ar 



r Jo V 4r r tj v ; 



r Jl + l/(Ar) 2 



r > r , A 1 . 



(A.3) 
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Now, depending on the relationship between A and ro, we obtain 

(A.4) 



^ In + A, r > r > A~ ] 



4 ln 4f' r>A- 1 >r , 



which gives Eq. (|3.10| ). 



APPENDIX B: ANOTHER REPRESENTATION FOR THE RETARDED 

GREEN'S FUNCTION 

Here, we obtain another representation for the retarded fermion Green's function that is 
more convenient for the derivation of the spectral density. Recall that when the imaginary 
part of G(w, k) is nonzero, \i + \/uj 2 — p 2 > and V\ > 0, v 2 < 0. This allows one to 
transform the analytically continued (by means of Eq. ( |4.18| )) integral 

i.r^, ,, 1 *;"'-' (B.D 

Jo [{u + v 1 ){u + v 2 )] a 



from Eq. ( |4.14|) as (for a < 1 



+ -U)] a J\V2\ [(«+ ui)(u - H)] Q wf 

T(l - (a, a, 1 - a; 1; « 2 ) + A 9^4 *1 f 1, a, 1 - a; 2 - a; - V (B.2) 
V V\ J \v2\r{2 — a) V v 2 v 2 J 



The first Appell function in (B.2) can be reduced to the hypergeometric function using the 



identity |4(| taht is valid for 7 = (3 + (3 1 

L-uT'VFi (a. 6: 6 + 3': 

i-y 



F 1 (a, /3, (3', (3 + (3'; x, y) = (1 - y)-« 2 Fi (a,/3;/3 + /3'; f— ^ J . (B.3) 



Thus one obtains 



■i)«r(a)r(i - a) / u 2 (l - Ul ) \ 1 r(l-a) 

-2-ri I a, a; 1; — - ) + - — rTT7 ^ x 



[ui(l-u 2 )] a V wi(l-w 2 )/ |« 2 |r(2-a) 

FxIW-a^-aAl); < 1, - < 0, - < 0. (B.4) 

M 2 M 2 / — «2j «2 U 2 



This completes the derivation of Eq. (|5.2j) . 



APPENDIX C: THE CALCULATION OF THE DENSITY OF STATES 



Introducing 



k 2 /2m ,11 H/ 

2/o = — — , (C.l) 



/i + y/uj 2 - p 2 ' 2m£ 2 _ /i + Vcj 2 - p 2 /i + Vu; 2 - p 2 
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and substituting ( |5.5|) in (|6.1|) , we can write 

| sgno; 6{u 2 — p 2 ) 



N(u) = N 

f 

Jo 



r(a) 



r(l — a) \mr 2 
dy 

2' 2 ' ' 



[{y + b- l) 2 + 46] Q / 2 



/a 1-a 

2-^1 77; ^ ' 



(^On^ + v/^-p 2 ) 1 -^ 

4y \ 



(y + b- l) 2 + 46 



x#(p + y^ 2 - p 2 ) - (^ 2 - p 2 -> -^ 2 - p 2 ) 

We now consider the integral from ( |U.2| ), 



4y 



/o [(?/ + 6-l) 2 + 46] Q / 2i " ± V2 2 
Using the relation f40| 



(y + 6- l) 2 + 46, 



2 Fx(a, 6; c; «) = (!- 2) a 2 Fi (a, c - 6; c; j^y) 



it can then be rewritten as 



dy 



(y + b + iy 



a 1 + a 

72-ri 77) — 77 — i i ) 



4y 



(2/ + 6+1) 2 



Replacing x 



6+1 
y + b + 1 



111 



we obtain 



x 



6 + 1 



2/0 + 6 + 1 



The integral (|C.5|) diverges at the lower limit as x — > or equivalently as y 
handle this we can write 



1= (b+l) l - a ^ 

J Xn 

= (6 + 1) 



dxx 



a-2 



J^;l;+H]_ 1 + 1 



6 + 1 



f 1 — T a ~ l 



a — 1 



To calculate the last integral in (|C.7|) , we rewrite it as 



E = lim / dxx 1 1 

y-*a-l Jo 



2Fl fei±^ ;1; Mi^)Ui 



2' 2 ' ' 6+1 y 
For 7 > we can compute the integral with the help of the formula (2.21.29) 



x a 1 (y — x) /3 1 2 -Fi(a, 6; c; ujx(y — x))dx 

o+fl-iD/ ox p / 7 o « + /?a + /5 + l ^y 4 ' 
y Q+/J l B(a,P) 4 F 3 [a,b,a,P;c,——, ; — — 



y, Rea, Re/? > 0, 
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so that 



E= lim ( B ( 7> i)^ s f« 1±H )7) 151,1 + 1, I±2 ; M_ I 
w ' ;4 3 V2' 2 ' ' ' 2 ' 2 '6 + 1/ 7 



S(a- 1,1)4^3 
iF (a - 1; 



2' 2 
a 1 + a 



a — 1 



2' 2 
1 



a — 1, 1; 1 



a 1 + a 1 



'2' 2 '6+17 a - 1 



1 



6 + 1/ a— 1 1 — en 



1-a 



6+1 



> 0. 



Thus, for the integral ( |(J.3[ ) we find 

1 



[(y + b+l) 1 - a -b 1 -" i 



1 — a 



Now substituting (|C.11|) into (|C.2|) we obtain 



N{u) = N 



a-1 



r(2 - a) \mrl 



sgnc^cu 2 - p 2 ) (A)n(/i + yjrf-p*) 1 -- 



(y + 6 + l) 1 "" - 6 1 "" + V^ 2 - P 2 ) - (V^ 2 " P 2 - " 



c^ 2 - p 2 ; 



Finally, replacing y and 6 in (|C.12| ) by expressions from ( p.l|) we arrive at Eq. ( 
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FIG. 1. Plot of the spectral function A(u, k) as a function of lo in units of the zero temperature 
gap A for k < k F , k = k F and k > k F at T = 0.99T B kt- 




FIG. 2. Plot of the spectral function A(lo, k) as a function of lo in units of the zero temperature 
gap A for k < k F , k = k F and k > k F at T = 1.043Tbkt- 




FIG. 3. Plot of the spectral function A(u>, k) as a function of lo in units of the zero temperature 
gap A for k < k F , k = k F and k > k F at T = 1.088T B kt- 
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FIG. 4. The density of states N(lj)/N at T = 0.99T BK t- 
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FIG. 5. The density of states N(u)/N at T = 1.043T B kt- 
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FIG. 6. The density of states N(u)/N at T = 1.088T B kt- 
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